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Abstract 

In a class of direct gauge mediation with a perturbatively stable SUSY breaking 
vacuum, gaugino masses vanish at the leading order of SUSY breaking F-term. We 
study the allowed parameter space of the gauge mediation models. By imposing a 
Tevatron bound on the lightest chargino mass m~± ^ 270 GeV and a warm dark- 

matter mass bound on the light gravitino mass m 3/ / 2 ~ 16 eV, we find that almost all 
the parameter space is excluded. Near future experiments may completely exclude, 
or possibly discover, the scenario. 



1 Introduction 



Low scale direct gauge mediation pQ is very attractive. It can achieve a very light gravitino 
mass m 3 /2 ^5 16 eV, satisfying the constraint from cosmology [2]. SUSY breaking vacuum 
can be sufficiently stable, even if not absolutely stable. Furthermore, low scale gauge 
mediation can be tested in future experiment at the LHC [3]. 

An explicit model of direct mediation with dynamical SUSY breaking was first con- 
structed in Ref. [I] (see also Ref. [5]) using the Izawa-Yanagida-Intriligator-Thomas (IYIT) 
model [6j as a SUSY breaking sector. Then, after the discovery of the Intriligator-Seiberg- 
Shih (ISS) metastable SUSY breaking model [7], many works have been done [8] to con- 
struct a direct mediation model by gauging a flavor symmetry of the ISS model. However, 
in those works it was observed [H El [8] that the Minimal Supersymmetric Standard Model 
(MSSM) gaugino masses are suppressed. That is, the gaugino masses vanish at the lead- 
ing order of a SUSY breaking F term, in an expansion in F/m 2 where m is a messenger 
mass scale. Later, the vanishing of the leading term in the gaugino masses is shown [9] 
in general set up. If a model can be described at low energies by some weakly coupled 
O'Raifeartaigh type effective field theory and the SUSY breaking vacuum of the model 
is the lowest energy state of the effective theory, the leading term of the gaugino masses 
vanish. 

To obtain the gaugino masses comparable to sfermion masses, one has to choose one 
of the following two possibilities. The first possibility is to give up the stability of the 
vacuum and go to a higher metastable vacuum [TO]. Such a vacuum becomes more 
unstable as we lower the gravitino mass (see Ref. [11] for a detailed study in the case 
of a so-called minimal gauge mediation). The other possibility is to have a messenger 
mass scale m and a SUSY breaking scale \fF to be comparable, i.e., m ~ yF. Then, 
higher order terms in the gaugino masses (of order 0(F 3 /m 5 )) are not so small compared 
with the sfermion masses (of order 0(F/m)). The condition m ~ \fF is also required to 
achieve the light gravitino mass. 

In this paper, we study calculable low scale direct gauge mediation in which a SUSY 

1 Note that this metastability is the one which is present even in the low energy effective field theory at 
perturbative level. The original ISS model is considered to be "stable" in this sense. (Non-perturbative 
dynamical SUSY restoration is not taken into account.) 
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breaking vacuum is stable at least at perturbative level. Then the gaugino masses are 
suppressed due to the vanishing of the leading order term. There are two important 
constraints on such a gauge mediation scenario; the constraint on the lightest chargino 
mass m-± ^ 270 GeV from the CDF collaboration at the Tevatron collider [12] (and also 
the recent work Ref. [13]), and the constraint on the light gravitino mass m.3/2 ^ 16 eV 
from the warm dark matter bound [5] a These constraints give lower and upper bounds on 
the SUSY breaking scale, respectively, so it is not obvious whether there is any parameter 
space consistent with the bounds. The aim of this paper is to study the allowed parameter 
space. We find that there is almost no parameter space consistent with the above bounds. 
There are some points in the parameter space that are near the margin of the bound, 
and at such points the gluino is quite light, nig ^ 800 GeV. Future experiments at the 
Tevatron (with enough integrated luminosity) or at the LHC may decide whether the 
scenario studied in this paper is relevant to nature or not. 

This paper is organized as follows. In section |2} we classify low scale direct gauge 
mediation models with a stable vacuum. Then in section [3] we calculate the maximum 
value of the lightest chargino mass m~± and the gluino mass m g under the condition 
that Yukawa interaction of the messenger sector is perturbative up to the GUT scale. 
We compare the result with the experimental bound. In section H] we consider a direct 
mediation in the ISS model, as an exceptional model in which the condition of section |3] 
is not satisfied. The last section is devoted to conclusions and discussion. 

2 Classification of direct mediation models 

In this section, we generalize the model discussed in Refs. [HE]. A direct gauge mediation 

model using the Intriligator-Seiberg-Shih (ISS) model will be discussed in section HJ We 

take the IYIT model [6] as a SUSY breaking sector for concreteness, but our result does 

not depend on this choice and any type of generalized O'Raifeartaigh models may work 

2 The gravitino mass bound is not applicable for a heavy enough gravitino mass and/or a low enough 
reheating temperature. But note that such a heavy gravitino mass (i.e. large SUSY breaking scale F) 
means heavy (i.e. well beyond 0(1 TeV)) sfermion masses due to the F/m 2 suppression of the gaugino 
masses. Furthermore, a low reheating temperature is inconsistent with many Baryogenesis scenarios, so 
we do not consider such a case in this paper. 
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as well. At low energies, the effective superpotential is given by 

W~A 2 Z. (1) 

Here, Z is a singlet chiral superfield, and A the SUSY breaking scale, respectively. We 
introduce Np flavors of messenger multiplets which transform as 5 and 5 under the GUT 
gauge group S77(5)gut- We will write the messenger quark multiplets as ^d,i,^d,i {i — 
1, • • • , Np) and the messenger lepton multiplets as 

We assume that the superpotential of the messenger sector is given as follows : 

W = A 2 Z + £ ^Jlf^Z)^ (2) 

X=d,l i,3 

(M x>ij {Z) = m Xtij + k Xtij Z) . 

We require that the model satisfies the following conditions. 

1. The SUSY breaking vacuum with = = is stable. 

2. The model has an R-symmetry. 

3. A D term of U(l)y is not generated by messenger loops at 1-loop level. 

4. There is no CP violation in the model. 

5. The Standard Model (SM) gauge couplings and the messenger sector Yukawa cou- 
plings are perturbatively small up to the GUT scale, 2 x 10 16 GeV. 

The first condition in particular requires that det(m + kZ) = detm [9]. In order for 
the messengers to be massive, we should impose that det(m + kZ) ^ 0. Then, detm ^ 0. 

The second condition is imposed to maintain the SUSY breaking [T4] (see also Ref. |15j). 
It is known that dynamical SUSY breaking models often have an R-symmetry. The IYIT 
model has an R-symmetry, and we assume that the messenger sector also respects that 
symmetry. Without the R-symmetry, introduction of "generic" Yukawa couplings between 
Z and the messengers restores SUSY [TJ]. However, note that this second condition does 
not completely ensure the first condition imposed above, because of the existence of a 
runaway behavior. See below. 

As to the third condition, if a non zero U(l)y D term is generated at 1-loop, the 
squared masses of the sfermions become negative. We introduce a messenger parity to 
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avoid such a danger [16] (see also Ref. p2]). The messenger parity is a symmetry that 
transforms the messenger chiral superfields and the U(l)y vector superfield V as 

follows § 



^ ^ £ UijVj, ^ ( 3 ) 

3 3 

where U and U are some unitary matrices. Owing to this symmetry, the dangerous D 
term is not generated by messenger loops. 

In order for the fourth condition to be satisfied, we assume that m x and k x are real 
for the time being. However, we will show that m x and k x can be taken real without loss 
of generality in the model which can make the gaugino masses maximum. 

The last condition requires that Np < 4, otherwise the perturbative gauge coupling 
unification is lost |18j . 

We investigate M(Z) satisfying those conditions. First, we simplify the messenger 
parity transformation by changing the basis of the messengers and redefining the trans- 
formation. We define ^[ = an d ^ = ^2j(XU)ij^fj, where X is some unitary 
matrix. Then, the messenger parity transformation is given by 

% % ^(XUUX^)^, V -V. (4) 

3 

Since UU is a unitary matrix, we can take X to diagonalize UU. Then the messenger 
parity transformation is given by 

% -> % % -> e^%, V -> -V. (5) 

When we operate the messenger parity transformation twice, ^' and \P' are transformed 

as, 

% -> f < ->■ e^tfj. (6) 

This transformation is also a symmetry transformation of the model. Then, in this basis, 

m . . = e m+^ mijj h . = e <W.+ft)^, (7) 

3 Of course SU(3)c and SU(2)l should also be transformed in an appropriate way. 
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We classify as follows : 

ft = 7T (i = l,...,^), (8) 
-7T<A<7T (i = JV w + l,...,jV F ). (9) 

For 1 < % < N n and + 1 < j < N F , < (3i + f3j < 2tt, then = rriji = kij = kji = 0. 
Therefore m and k are block-diagonal matrices which are given by 

For JV* + 1 < i, j < N F , -2ir < fa + < 2vr, then, & + = or = k {j = 0. Now, 
we show that the following transformation is a symmetry transformation of the model. 

%^Y,% %^Y,% (11) 

where 

y _ \ 1 i = l,-",iV w , 2 x 

The proof is as follows. Under this transformation, 

E E k-//» ;j -v>W- (is) 

Because + — for i, j — N n + 1, • • • , JVp for which m^- 7^ 0, 

™ = ( m * row e-«^>/> ) = ( m * my ) • < 14 > 

Therefore, the mass term is invariant. It is easy to check that the Yukawa interaction 
term is also invariant. 

By using the Y transformation, we define a new messenger parity transformation. 
After operating the original messenger parity transformation, we operate the Y transfor- 
mation. Under this transformation, ty', and V are transformed as 

(< = 1,...,JV W ), 
%^e~ iM2 %, %^e^l 2 % {i = N n + l,---,N F ), 

V ->■ -V. (15) 
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From now on, we call this new transformation as the messenger parity transformation. 
By using appropriate phase rotation, we can take a new basis and in which the 
messenger parity transformation is given as follows : 



Vi^^i (i = N n + !,•••, N F ), 



V -> -V. 



(16) 



Because m and k are block-diagonal, we consider the cases N n = and N n = Np sepa- 
rately. Then, we have U — 1 and U = ±1. In more general cases, the messenger sector is 
a direct sum of these two cases. 

Let us next consider about R-symmetry. We assumed that there is an R-symmetry. By 
using this R-symmetry and the messenger parity, we can define the following symmetry 
transformation, 



where R is the generator of the R-symmetry on \l> and R is on ty. We have used U = 1 and 
U — ±1. Under this transformation, the superpotential W has R-charge 2. Therefore this 
transformation also generates an R-symmetry. Then, we can define a new R-symmetry 
transformation as the sum of the original and the above R-symmetries, 



It is obvious that this new R-symmetry commutes with the messenger parity. There- 



symmetries. From now on, we call this new R-symmetry as the R-symmetry of the model. 

We assign R-charges as -R(^i) < R($ 2 ) <■• < R($n f )- We now show that R($i) + 
R(V Np -i + i) = 2. If R{^i) + R(y NF _ i+ i) < 2, • • • , *j cannot have mass terms with 

• • • , ^N F -i+i- Then the mass matrix is of the form, 




(17) 
(18) 



q, _^ e ia{R+R)/2q,^ q, _^ e ia(R+R)/2^_ 



(19) 



fore, we can assign R-charges as R(^i) = R(^i) by simultaneously diagonalizing these 




(20) 
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where Oi t N F -i+i is a i x (Np — i + matrix with all components equal to 0. This matrix 
has detm = 0. One can also show that if R(^i) + R(^N F -i+i) > 2, then detm = 
0. This contradicts with the stability condition discussed above. Thus, we must have 
R(Vi) + R(^ NF . i+1 ) = 2. 

We write R(^i) as a. If W includes an interaction term kZ^!{^> jv>-i+ij we have 
R($?N F -i+i) = _ a , and R($i) = a+2. Then if W includes an interaction term kZ^i^f^p-j+i, 
R(^N F -j+i) = —a — 2 and R(^j) = a + 4. In this way, there are ^ which have R-charges 
a, a + 2, • • •, and * which have R-charges —a + 2, —a, • • •. We call them as a-group. If 
there are some \I/ which are not included in a-group, we write the smallest R-charge of 
them as 6. Such \1> make up o-group. Therefore, all * are classified into the groups, 
a-group, 6-group, c-group and so on. The components of M X (Z) are nonzero only within 
each group. Note that if we take one group, R(fy i+1 ) — R(^i) = or 2. 

Some groups have a messenger parity by itself. If not, a group (assume it to be a-group) 
has an "anti-group" (which is (— a — 2n)-group, where (— a — 2n) is the smallest R-charge 
of \& contained in a-group) and they make up a pair to exhibit the messenger parity 
symmetry. However, we can combine these two groups into a single group by defining a 
new R-symmetry. Since the components of M X (Z) are nonzero only within each group, 
there is a U(l) symmetry which assigns charge x to ^ in a-group, —x to ^ in a-group, y 
to ^ in (—a — 2n)-group, — y to ^ in (—a — 2n)-group, and to other fields, where x and 
y are arbitrary numbers. Then, by using this U(l) symmetry with x = —y = —(a + n), 
we can define a new R-symmetry which assigns charge — n, — n + 2, •••,n + 2to\I / and 
n+2, n, • • • , — n to ^ in both of the groups. Thus, we can combine a-group and (— a+2n)- 
group into a single (— n)-group, which has the messenger parity by itself. Therefore, we 
can assume that the messenger sector is a direct sum of groups with each group having 
the messenger parity by itself. Note that we can maintain i?(\l/j + i) — R(^i) = or 2 in 
the combined group after appropriately rearranging \l/ such that R(&i) < R(^i + i). 

In summary, we have obtained R(^i) = R(^i), R(^i) + R(^N F -i+i) = 2, and 
R(^f i+ i) — R(^i) = or 2 in one group. If we determine an R-charge assignment of 
the messengers following the above conditions, we can get the form of M X (Z). 

First, we consider the case U — U = 1, and there is only one group. As noted above, 
more general cases can be expressed as a direct sum. The messenger parity transformation 
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is given by % <H> and thus M X (Z) = M%(Z). 

• Np = 1 model : 

In this model, M X (Z) = m, that is, the messenger has no Yukawa interactions. 

• Np = 2 model : 
(R(V 1 ),R(V 2 )) = (a,a + 2), a = 0. 

(kZ m \ 
m J' 

• N F = 3 model : 

i2(* 2 ), #(* 3 )) = (a, a + 2, a + 4), a = -1. 

kZ m! 
\ m 

• Np — 4 model A : 

#(* 2 ), i?(* 3 ), ^(^4)) = (a, a + 2, a + 4, a + 6), a = 

/ /cZ m \ 

k'Z m! 
kZ m! 
\ m ) 

N F = 4 model B : 

i2(* 2 ), i?(* 3 ), i?(*4)) = (a, a + 2, a + 2, a + 4), a = -1. 

/ fciZ fc 2 Z mi \ 

k\Z m 3 m 2 
A; 2 Z m 2 m 4 
\ mi / 

N F = 4 model C : 

i2(* 2 ), 3), i?(*4)) = (a, a, a + 2, a + 2), a = 0. 

/ fc 2 Z /ciZ m 3 mi \ 

fciZ k 3 Z m 2 m 4 

m 3 m 2 

\ mi m 4 / 



(21) 



(22) 



(23) 



(24) 



(25) 
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In fact, some of the above models have a runaway direction in the tree level potential. 



See |Appendix A| for details. We must require that km 1 k = for the model to have 



no runaway direction. Np = 3 model has a runaway direction. Np = 4 model A has a 
runaway direction unless k' = or k = 0. By unitary rotations, we can take k 2 = in 
Np = 4 model B. Then this model has a runaway direction unless m 4 = 0. Np = 4 model 
A is included in Np = 4 model B (if k' = 0) or Np = 2 model (if k — 0). There is no 
runaway direction in Np = 4 model C. 

So far, we have considered the cases in which there is only one group. We can consider 
a direct sum of two N F = 2 models, but this is included in Np = 4 model C. 

Next, we consider the case U = 1, U — — 1. The messenger parity transformation is 
given by ^ -> ^ -> and thus M X (Z) = -M%(Z). In this case, Ajr is even, 

otherwise detM x (Z) = 0. 

• Np = 2 model : 

(i2(*a),i2(* 2 )) = (a,a + 2), a = 0. 



-m 



(26) 



In this model, the messengers have no Yukawa interactions. 
Np = 4 model A' : 

(R(Vi), R{®2), R(V 3 ), R(V*)) = (a, a + 2, a + 4, a + 6), a = 



-2. 



/ 

— kZ —ml 
\ —m 



kZ m \ 



m 



/ 



(27) 



N F = 4 model B' : 

(R(Vi), R(* 2 ), R(y 3 ), R(^4}) = (a, a + 2, a + 2, a + 4), a 



/ fciZ k 2 Z mi \ 

—k\Z m-i 

-k 2 Z -m 2 

V -mi / 



(28) 
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By unitary rotations, we can set k\ = 0, i.e. 



/ 

—k 2 Z —m 2 
\ -mi 



k 2 Z mi \ 
m 2 

J 



Np = 4 model C : 

i2(* 2 ), #(^s), i?(* 4 )) = {a,a,a + 2,a + 2), a = 0. 



( kZ m 3 mi \ 

— /cZ m 2 m 4 

— m 3 — m 2 

\ —mi — m 4 J 



By unitary rotations, we can set mi = m 2 = 0, i.e. 

/ kZ m 3 \ 

—kZ m 4 
-m 3 

— m 4 



(29) 



(30) 



(31) 



All these models are in fact equivalent, and are included in Np = 4 model A with k' = 0. 
After all, models which satisfy all the conditions are the following three models. 
N F = 2 model : 



kZ m 
m 



(32) 



Np = 4 model B 



Np = 4 model C : 



( kZ m\ \ 

kZ m 3 m 2 
m 2 

\ mi / 



(33) 



/ kZ m 3 mi \ 

k'Z m 2 m 4 
m 3 m 2 
\ mi m 4 / 



(34) 
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where we diagonalized the Yukawa couplings by unitary transformation. 

Let us comment on CP violation. The parameters of Np = 2 model and Np = 4 model 
B can be taken real without loss of generality. Then, there is no CP violation. However, 
Np = 4 model C has CP violation in general. To ensure the absence of CP violation, we 
assume that all parameters in the model are real. 

3 Upper bound on gaugino masses 

To avoid a cosmological problem, we want to achieve a light gravitino mass m 3 / 2 ^ 16 eV. 
But the lightness of the gravitino mass means smallness of F, which is the F component 
of the chiral supermultiplet Z. Then the gauginos become light. 

To get large gaugino masses, it is important to have a large Yukawa coupling k, by 
the following reason. The messenger mass spectrum depends on the combination kF, and 
the gaugino masses can be roughly written as 



where f(kF/m 2 ) is a dimensionless function of kF/m 2 . In gauge mediation, f(kF/m 2 ) 
is bounded from above 3, f(kF/m 2 ) < / max . Then, the gaugino masses are also bounded 
from above as 



where Mpi ~ 2.4 x 10 18 GeV is the reduced Planck mass, and we have assumed that the 
F term of Z is responsible for all the vacuum energy. Therefore, we want large Yukawa 
coupling constants to achieve large enough gaugino masses. However, we can not take 
arbitrary large values for the Yukawa couplings, because the Yukawa interactions are 
not asymptotic free. Then, we have to know the allowed parameter region in which the 
Yukawa interactions are perturbative up to the GUT scale. 

Such a parameter region is determined by renormalization group (RG) equations. The 
RG equations of the Yukawa couplings in each model of the previous section are given 

4 Generically, the maximum value is achieved when masses of some scalar components of messenger 
chiral fields become almost zero. 




(35) 




(36) 
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Appendix B By using the anomalous dimension of Z and the messengers, the RG 



equations of the Yukawa couplings are given by 

dk k 

-^=1^(7Z+7« X +7* X )(X = 40, (37) 

where t = log(/i/yU ), with fi the RG scale. 7^ x and 7^ depend on fc 2 for each x — d, I, 
but is proportional to 3k% + 2kf. This means that the RG equation of k\ depends on 
kd through jz- Then, if we make kd smaller, we can get larger k%. 

First, we investigate the maximum value for the wino mass by setting kd = 0, and see 
whether that value can exceed the CDF bound [121 13] Tn^± ~ ^70 GeV (we will discuss 
on this bound later). If kd 7^ 0, the wino mass bound becomes smaller than in the case 
k d = 0. In the case k d = 0, the RG equations in each model are given as follows : 

Np = 2 model : 

dh - h Ukf-3g!- 3 -gl). (38) 



dt 16vr 2 V 1 1 ^ 



Np = 4 model C: 



4Af + 2Af - 3gl - ^ , (39) 



= ( A i3 , nu'2 Q „2 3 ^2 

dt 167T 

dk'i k[ f 2 /2 2 3 



a i6.n 2fc ' + <- 3 *-5^- < 40 > 

The RG equation for Np = 4 model B is obtained by setting k[ = ki in Eq. ( 139]) . 

By using these RG equations, we can determine the upper bounds on k\ (and k[). We 
require that the Yukawa coupling constants at the GUT scale are less than Att, but our 
result does not strongly depend on this value as long as it is large enough. 

The upper bounds on the Yukawa couplings at the messenger mass scale in each model 
are as follows. 

• Np = 2 model : 1.06 

• Np = 4 model C : Figure [1] 

In these parameter regions, we calculate the wino mass. We do not specify the mechanism 
generating the R-symmetry breaking vev (Z) (see e.g. [IE]), and treat ki, mi and (Z) as 



free parameters. In the calculation, we used the formulae in Appendix C By using the 
upper bound on ki, we can calculate the upper bound on the wino mass. 
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Figure 1: The upper bound on the Yukawa coupling constants for Np = 4 model C 
obtained by solving the RG equations ( l39lT4T)i) . The bound for Np = 4 model B can be 
obtained by setting k[ = ki. 



The result for the wino mass using «2 = Oiiimz) is, (we will give physical pole masses 



later) 



2 model : m w < x 0.24 ~ 160 GeV x ( 



V16 eV 



4 model : m w < — \[¥ x 0.43 ~ 300 GeV x ( 



47T 



V16 eV 



1/2 



1/2 



where we have taken into account the upper bound on the gravitino mass 7113/2 ^ 16 eV. 
The upper bound is the same for both Np = 4 model B and C. 

The actual upper bound on the wino mass (i.e. when kd ^ 0) is lower than those 
results. Therefore Np =2 model is clearly in conflict with the CDF bound, while Np = 4 
model is marginal. Let us focus on this case. 

In Np = 4 model, we found that the gaugino masses can be maximized when M X (Z) 
is given by 



M X (Z) 



( kyZ 



k x Z m x 

777.V 



m x \ 



(41) 
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Figure 2: The upper bound on the Yukawa coupling constants kd and k\ obtained by 
solving the RG equations (142f43p . 



In the rest of this section, We assume that Np = 4 and M x is of this form. Fortunately, 
in this model, m x and k x can be taken real without loss of generality, so there is no CP 
violation B 

Now we have six parameters, kd, fcj, m^, mi, (Z) and 777.3/2- kd and ki are restricted by 
the condition that the Yukawa interactions are perturbative up to the GUT scale. The 
RG equations of the model are as follows : 



dk d k d ( 2 2 16, 4 2 



8kj + Akf - ^gi - i-gi) , (42) 



dt 16vr 2 V a 1 3 " 6 15 

dkl h f c 2 , c l2 q„2 3 „2 



* 16^V 6 « + W '- 3 *-5^' (43) 
First, we calculate the allowed parameter region of k d and k\. The result is shown in 
Figure |2j By using it, we calculate the sparticle masses in whole the parameter space of 
kd, h, rrid, mi and (Z) with 7773/2 < 16 eV. For the pole mass calculation, we used the 
program S0FTSUSY 2.0.18 [20] We assumed that sign(/x) = +1 and tan/3 = 10, where 



5 We thank T. T. Yanagida for pointing this to us. 

6 We thank S. Shirai for preparing and operating the program. 
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Figure 3: The upper bound on the lightest chargino mass m~± and the gluino mass ms. 
The CDF bound m-± ^ 270 GeV is also shown. 



/x and tan/3 are the usual MSSM parameters (see e.g. Ref [21]). but the result for the 
lightest chargino and the gluino masses is nearly independent of this choice. We show the 
upper bound on the lightest chargino and the gluino masses in Figure [3J 

Let us discuss the experimental bound on the models discussed above. In the CDF 
collaboration at the Tevatron collider [12], a bound on the soft masses was obtained 
by assuming a minimal gauge mediation model. Their parameter choice is : N m = 
l,M m = 2A m ,tan/3 = 15, and sign(/i) = +1, where N m is the messenger number, M m 
the messenger mass scale, A m the scale by which the gaugino masses are given by = 
N m (ai/4:7i)A m . They obtained a bound m^o > 149 GeV for the lightest neutralino Xi 
(which is the next to lightest supersymmetric particle (NLSP)), and the bound can be 
translated into a bound on the mass of the lightest chargino m -± > 290 GeV in that model. 
In these parameters, Xi is bino-like, and the next-to-lightest neutralino xH anc ^ ^ ne lightest 
chargino Xi are wino-like. Dominant production processes are q + q — > (7/Z )* — > X1X1 
and q + q — > (W ± )* — > X1X21 an d they search for a di-photon+missing energy event which 
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comes from the decay Xi ~~ jG (G is the gravitino) at the end of the cascade decays of 
produced particles. 

In the models considered in this section, Xi is the NLSP, and xt > X2 an d Xi are wino- 
like and bino-like respectively, except for a special parameter space discussed below. This 
is because the gauginos are somehow lighter than other sparticles due to the suppression 
of the gaugino masses discussed in the Introduction. In Ref. [13], the CDF data was 
reanalyzed in the context of a general neutralino NLSP. They obtained a conservative 
bound on the lightest chargino mass m-± ^ 270 GeV in the case m_g < <C //. Note 
that the production cross section depends on the mass of xt and X21 which is almost equal 
to the wino mass mw in the present model. The bound m~± ^ 270 GeV exclude almost 
all the parameter space, as we can see from Figure EJ Furthermore, the Tevatron may 
achieve a stronger bound [13] m-± ^ 300 GeV with 10 fb _1 of integrated luminosity. Even 
if not excluded completely at the Tevatron, one can see from Figure [3] that the gluino 
mass should be quite small, m~ g < 800 GeV, to be discovered at the LHC. 

Let us briefly discuss other possibilities for the NLSP which is not a bino-like neu- 
tralino. If the NLSP is not a bino-like neutralino, the branching fraction to photons is 
suppressed compared with the bino-like neutralino case. One candidate for the NLSP is 
the stau, f\. We checked that the stau cannot be the NLSP for tan/3 as large as 50. A 
wino-like or higgsino-like neutralino NLSP is possible in the present model (e.g. if k\ <C kd 
or ki 3> kd, respectively), but then the mass is small, m^o ^$ 90 GeV. This case is also 
excluded by the CDF bound [13]. Also there is a region where the lightest chargino is 
the NLSP, but the mass is small, m^± ^ 90 GeV. Thus this case is excluded by LEP 
bounds [22J. Other sparticles are too heavy to be the NLSP. 

Finally we comment about possible effects on the RG equations from interactions other 
than the messenger sector Yukawa couplings. If Z has interactions with other fields in the 
SUSY breaking sector, as in the IYIT model, 72 presumably receives positive contribution 
as long as Z is a gauge singlet. This makes the messenger Yukawa couplings smaller than 
that obtained above, leading to a further suppression of the sparticle masses. Introducing 
new singlets which couple to the messengers also make the Yukawa couplings smaller. 
On the other hand, the model Eq. ( HT|) has an SU(2) global symmetry, which can be 
gauged. The gauge interactions give negative contribution to the messenger anomalous 
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dimensions, and this makes the Yukawa couplings larger. Unfortunately, we found that 
the effect does not drastically improve the situation, maximally 12% increase of the wino 
mass. It is a model building challenge to find a mechanism which can make the Yukawa 
couplings larger. One such mechanism is discussed in the next section. 

4 Direct mediation in the ISS model 

In the previous section we have studied direct gauge mediation models under the condition 
that the Yukawa couplings do not blow up below the GUT scale. Let us next consider a 
direct gauge mediation in the ISS model [7]. In this case, Yukawa couplings are generated 
dynamically. Then we can take the "UV cutoff scale" of the Yukawa couplings, A cut , 
to be the dynamical scale of the ISS model. The Yukawa couplings should not blow up 
only up to the scale A cut . As we lower the cutoff scale A cut , the Yukawa couplings at 
the messenger scale can become larger, leading to larger gaugino masses. Note that this 
dynamical generation of the Yukawa couplings is an advantage of the ISS model over 
simply retrofitted [23] O'Raifeartaigh type models. 

The model we consider is based on an S'?7(5)hid massive SQCD with 7 flavors of 
quarks and anti-quarks at high energies Q. After the confinement of »S?7(5)hid, the low- 
energy theory is given as follows [21]. The matter content of the model is; mesons $j 
and (anti-)quarks (pf, (p x a . Here I, J — 1, • • ■ , 7 are flavor indices (part of which we will 
identify as SU(5)gut indices) and a = 1,2 a gauge index of the dual magnetic gauge 
group 5 , f/(2) mag . They are decomposed as 



Here we have decomposed the indices I (J) into / = {p, a} (J = {q, /?}), where p = 1, 2 
is a flavor index and we identify a = 1, ■ • • , 5 as the SU (5) gut index. 
At the scale A cut , the superpotential is given by 



7 In choosing this model, we have taken into account the Landau pole problem of the SM gauge 
coupling constants. See the discussion at the end of this section. 




(44) 



W 



ktr(ip^0) - tr(yu 2 $) 

k [tr(^Z^) + tr(x^'^ + X^') + ^(x Y X)} ~ tr(/z 2 $) 



(45) 
(46) 
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where we have assumed that there is a single Yukawa coupling constant k at the scale 
A cu t- Without loss of generality, we can assume that all parameters are real and positive. 
Because of the renormalization group effect, there are several Yukawa coupling constants 
at low energies, 



ktr(^Z^) -> tr 



( 



k' dl Z d + k dl l 3 ■ z d /Vs k dll Z dl \f V a 



(47) 



where Z d is an SU(3)c adjoint field, Z\ an SU(2)l adjoint field, Z d i, Z[ d bifundamental 
fields of SU{3) C x SU(2) L , and Z d , Z l singlets. Also, 

hr(xrt + xW') -> *4tr(x$i$j + x*^) I ^Mx*!*! + X*i*U (48) 
kti( X Yx) ->> A; 3 tr( X Yx). (49) 

The RG equations for these couplings are given in Appendix B The term tr(/i 2 $) is, in 



the above notation, given by 

tr(/7 2 $) = ii\ tr(Y) + V3fi 2 Zd Z d + v^/ii^. (50) 
At a SUSY breaking vacuum, x, x develop vevs, 

(x;) = <X> *£, (XS> = (X) SI, (x) = ^ V (51) 
Then, we obtain a gauge mediation model 

W = £ [m x tr(tf' x # x + * x *' x )] + k dl Z d tr(tf d tt d )/V3 + fc,^ tr(tf ,)/>/2 

-V3//| d Z d -V2//|^ + ---, (52) 

where m; = (x) an d m d = k d 2 (x)- This is similar to the model considered in the 
previous section. The vacuum energy (in global SUSY) is given by V = 3/if d + 2fi zl , and 
the gravitino mass by m 3 / 2 = W /y/3Mpi. We require that 777,3/2 < 16 eV. 

There is a subtlety in choosing the values of (x = d, I). If we impose that the 
messenger fields *ff, ^ , and \P' do not have tachyonic masses, /i| must satisfy 

ii\ x < k-\m\ = ^fy (x = d, I.) (53) 
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On the other hand, for the above vacuum to be a global minimum of the potential at tree 
level, it is necessary that ^z x satisfy 

/4 < My, (54) 



'2k 2 + k' 2 x 1/2 



^ < i d i; 2 di (55) 

V (Z1 / 



These conditions can be obtained by comparing the vacuum energy at \ ^ 0, \&j = 
0, = (the above vacuum) to the vacuum energy at x = 0, ^ ^ 0, ^ = or 



at x — 0, \l/z = 0, 7^ 0. See |Appendix D| for details. The RG equations make 
Eqs. f l54|55|) stronger than Eq. (1531) . Requiring Eqs. f )54ll55|) leads to a suppression of the 
gaugino masses. 

We calculate the maximum values of the lightest chargino and the gluino masses as 
a function of the cutoff scale A cut . We assume that all the Yukawa couplings and the 
magnetic gauge coupling g m are 47r at the scale A cut . We do not assume model building 
details and simply take the R-symmetry breaking vevs (Zj) and (Zi) as free parameters. 
See Ref. [8] for mechanisms generating such vevs. The result is shown in Figure H] and 
[5j We show the upper bound under either the condition Eqs. f)54|55p (solid line) or 
Eq. (loUj) (dashed line). However, for our purpose of investigating a stable SUSY breaking 
vacuum at perturbative level, only the solid line should be taken seriously. As discussed in 
section |3l the lightest chargino mass is bounded as m-± ^ 270 GeV from the CDF bound, 
which is possible only for a very small cutoff scale A cut /M mess ^ O(10). Since the ISS 
model is metastable, A cut /M mcss is bounded from below to suppress the vacuum tunneling 
rate of the model. Detailed numerical study is required to determine how we can lower the 
cutoff scale A cut . (See Ref. [7] for a crude estimate of the vacuum tunneling rate.) Note 
also that the gluino mass is small, nig 1z 800 GeV, for the cutoff scale A cut /M mess ^ 10. 

Finally, let us comment on the perturbative unification of the SM gauge couplings 
at the GUT scale. Direct mediation models in the ISS model have many fields charged 
under the SM gauge group. Then, the SM gauge couplings hit Landau poles below the 
GUT scale (or at least the precise perturbative unification is lost [18]) if we consider a low 
scale direct mediation with m 3 / 2 < 16 eV. Recently, it was shown that the perturbative 
unification can be maintained by appropriately modifying the theory [22]. In particular, 
a modified version of the above model is shown to be consistent with the perturbative 



20 



350 



300 



250 

> 

a 200 



150 



100 



50 



10° 10 1 10 2 10 3 10 4 10 5 10 6 10 7 10 8 10 9 10 10 10 11 



A cu/ M r 



Figure 4: The upper bound on the lightest chargino mass m-± in the ISS model as a 
function of A cut /M mess , where A cut is the cutoff scale and M mess the messenger mass scale. 
The solid line represents the upper bound under Eqs. (I54ll55p . while the dashed line under 
Eq. (|53l) . The CDF bound m ~± ^ 270 GeV is also shown. The vev (Zi) is chosen to 
maximize the lightest chargino mass. 
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Figure 5: The upper bound on the gluino mass rrig in the ISS model as a function of 
A cut /M mess . The solid line represents the upper bound under Eqs. ( 153|55|) . while the 
dashed line under Eq. (I53p . The vev (Za) is chosen to maximize the gluino mass. 
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unification [25]. In the modified model, the low energy effective superpotential of the 
model is slightly different from the one described above. It is interesting to study further 
in that case. See also Ref. [26J for another approach to the GUT unification. 

5 Conclusions and discussion 

In this paper we have studied the parameter space of low scale direct gauge mediation 
models with a (perturbatively) stable SUSY breaking vacuum. We have found that the 
CDF bound on the lightest chargino mass m~± ^ 270 GeV and the warm dark matter 
bound on the gravitino mass 777,3/2 ^5 16 eV severely restrict the allowed parameter space 
of the low scale direct mediation models, and almost all the parameter space is excluded. 
Even if the models somehow manage to satisfy the bound, the gluino mass is quite small, 
m g ^5 800 GeV. Such a light colored particle is supposed to be discovered at the LHC 
experiment. Furthermore, in the models studied in this paper, the NLSP is a bino-like 
neutralino, which decays to a gravitino and a photon \\ — >■ G + 7. Thus there is a di- 
photon + missing energy signal in a SUSY event, which makes the discovery more easy. In 
addition to the LHC experiment, future cosmic microwave background surveys may give a 
stronger upper bound on the gravitino mass [27], m 3 / 2 ^ 3 eV. Thus, the scenario studied 
in this paper should be completely excluded, or discovered, in near future experiments. 
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Appendix A Runaway direction 



We will investigate the condition that a direct mediation model has no runaway direction. 
The superpotential of the model is given by 

W = XA 2 Z + *(m + fcZ)#, (A.l) 

where det(m + kZ) = detm ^ 0. We assume that the Kahler potential is canonical. 
Then, the scalar potential is given by 

V = \XA 2 + ^k^\ 2 + \(m + kZ)^\ 2 + \^(m + kZ)\ 2 . (A.2) 

We define k and k as 

K = (m + kZ)^, (A.3) 

k = ^{m + kZ). (A.4) 

Because det(m + kZ) = detm ^ 0, (m + kZ) is invertible for arbitrary Z. Then, by using 
k and k, the scalar potential is written as 

V = \\A 2 + k(m + kZ)- 1 k(m + kZ)- 1 K\ 2 + \k\ 2 + \k\ 2 . (A.5) 

If the model has a runaway direction, in that direction each term must be arbitrarily 
small, so the norm of k and R must be arbitrarily small. In that direction, the first term 
in Eq. (1A.5I) can become small if and only if (m + kZ)~ 1 k(m + kZ)~ l depends on Z. 
Therefore, the condition that the model has no runaway direction leads to 

(m + kZ)~ 1 k{m + kZ)^ 1 = m~ l km~ l 

k = (m + kZ)m~ 1 km~ 1 (m + kZ) 

= k + 2km~ 1 kZ + km- 1 km' 1 kZ 2 . (A.6) 

Then, the model has no runaway direction if and only if km~ l k = 0. 

Appendix B RG equations 

The RG equations of the Yukawa coupling constants are given as follows. 72 is the 
anomalous dimension of Z. 



23 



N F = 2 model : 

dk d ( 1 / a 16 2 4 2 
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dk d ( 1 / 2 16 2 4 2 
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16vr 2 
fc 3 

167T 2 

16vr 2 



_t 2 + rt' 



3 T~> 2 
2' 



on 'dll 



+ 3*4 + 64 + 2A; 3 2 - 3g\ - ~g\ ~ 3 



Qk 2 d2 + 4k 2 2 + Qkl-3gl 



24 



where g m is the gauge coupling constant of the magnetic gauge group S77(2) mag . 

Appendix C Gaugino and sfermion mass formulae 



In this appendix, we will give formulae for the gaugino and sfermion masses □ at the 
leading order of the SM gauge couplings, when the low energy effective superpotential is 
given by 

w = E E (mffv^vW + kffz^vf) . (C.l) 

X=d,l i,j 

After integrating out auxiliary fields, mass terms of the messengers are given by 
J2 jV x) m /(x V (x) + h.c. 

h V^ ^ ) \ -k<x*(F z ) m'Wm'Wt I \ 0(x)t 



--mass 

X =d,l 



(C.2) 



where m'^ x > = mpv + (Z), ip^ and ip^ denote the messenger fermions, and <ffiu 
and <fi^ denote the messenger scalars. We can diagonalize the mass matrices by unitary 
transformations. There are unitary matrices U^ x \ and which diagonalize the 
mass matrices as follows : 

[/(x)t m '(xV(x) = M (x) ; ( C3 ) 

{x)t / m'Wt m '(x) - k W(F z y \ {x) _ (x)2 
H [-k<x)( Fz ) m '(x) m '(x)t ) R ~ M s , ^-4) 

where = diag(mj^, m^, • • •) and M^ 2 = diag(m^ x 1 ' )2 , m^ 2 , • • •)■ We can define 

mass eigenstates -0^ x / and as follows : 



(x) = R Wi ( f x) \ 



(C.5) 



^(x) = y(x)t^,(x) ; ( C6 ) 
^(x)t = [/ (x)t^(x)t. (c.7) 



For general mass formulae of weakly coupled gauge mediation, see Ref. 
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We define the upper half of i?W as A&) and the lower half of R<x> as fiW. and 
i?^ x ^ are n x 2n matrices, 

— (£)■ 

By using it, we get the formulae for the sfermion masses 



(C.8) 



and the gaugino masses 



(0 



m 9i = ^ ^1/2 + 5 m i/2 



^§3 



«2 (0 

27T mi / 2 ' 
«3 (d) 

27 m V2- 



(C.9) 

(CIO) 

(C.11) 
(C.12) 



In Eq. (IC.9p . C 2)S and C 3iS are the quadratic Casimir invariants, and Y s denotes the U(l)y 
hypercharge. Here uIq and rrt^X are given by 



rn 



(X)2 



n (v)2i fv")2 (y!2 , 

2 mfl log m£/ + m\y log m 



,(X)2 



(X)2 



,(X)2 
L s,i 



(X)2' 



III 



(X)2 



2 • • ' 



ml i Li 2 1 



Hi 



(X)2' 
'a,3 



in 



(X)2 
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1/2 



^(i W V (x) ) 13 (f/ Wt B (x) ) 



(x)2 (y) 
m s,i m f,j 



log 



m 



Wa- 



rn 



(X)2 
/j 



where Li 2 (x) is the dilog function 



oo 



fc=l 



i rft log(l-xt) 
£ 



(C.13) 
(C.14) 

(C.15) 



Appendix D Potential minimum of the ISS model 



Due to the RG effects on the Yukawa couplings in Eqs. (j47|48|49p . the potential of the 
ISS model is complicated. The potential from F-term equations is given by 



V 



k dl tr(tf d tf d ) - 3/4 d 



+ k 



1 2 
dl 
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+ 2 



1 _ 2 
- -l a tr(*,*,) 



+ 



+k 2 
+ •-• 



+ kj 2 V dX + |x* 



t.2 



(D.l) 



where dots denote terms which can be set to by appropriately choosing the meson vev 
($) (e.g. (<&) = 0). Here ^ is a 3 x 2 matrix, VP^ is a 2 x 3 matrix, and others are 2x2 
matrices. For simplicity, let us pretend as if k^, k\i and kdii are infinitely large. Then, 
the last line of Eq. ( ID.lj) must vanish, and we should consider the minimization of the 
first three lines under that condition. 

As a preparation, let us consider the minimization problem of a potential of the fol- 
lowing form, 



1 jfctr(^) -NfJ 2 2 + k' 2 



W> ~ ^lNtr(^) 



(D.2) 



where ip is a N x n' matrix and ip is a n' x N matrix with N > n'. We minimize Eq. (ID.2[) 
under the condition that rank('0^') = n with n < n' . By SU(N) unitary transformation, 
we can assume without loss of generality that 



ipifj 



A 
* 



(D.3) 



where A is a n x n matrix. It is easy to see that at a potential minimum, the * in Eq. (ID.3[) 
is zero. Then, 



V > — 

~ N 

= k' 2 



1 ktr(A) -Nfi 2 \k' 2 



A--l n tr(A) 



+ k' 2 



-^lN-ntr(A) 



1 2 

A l n tr(A) 

n 



tr(A)| 2 - 2fc/i 2 Retr(A) + Nfi 4 



(DA) 



From this form, we can see that the minimum is 



Vrnir 



N(N - n)k' 2 



nk 2 + (N — n)k 



72^ 



(D.5) 
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at 



.4 



n \ N Nn 



fc/^ 2 l n . 



(D.6) 



By using Eq. flD.5j) . we can see that the minimum value of the potential fID.ip is given 
by one of the following possibilities. 



V m 



3/4d + 2 /4 

3/4d + 2 /4 



Zl 



K ii +K n 



at rank(xx, Wj, W c 



S^^i 4 I o 4 , 4 
6fc^ 4 . 2fe^ 4 „ 4 



(2,0,0) 
(0,2,0) 
(0,0,2) 

(0,1,1) 
(1,0,1) 
(1,1,0) 



(D.7) 



where the rank of the matrices ^d^d) is determined by requiring that the last 

line of Eq. (ID.lj) should vanish. By solving the RG equations presented in Appendix B 
numerically, we checked that k d x < k' dl and kn < k[. Then, requiring that the first of 
Eq. ( ID. 71) is the smallest, we obtain 



2^2 



Vzd < 



2k dl + k' 



1/2 



(/1 



3k 2 



(D.8) 
(D.9) 



en 



Under these conditions, the vacuum with rank(xx, ^1^1, ^d^d) = (2, 0, 0) is stable. 

So far we have assumed that k d 2, ki 2 and k d n are infinitely large. It is obvious that 
Eqs. ( ID.8fD.9j) are necessary for the above vacuum to be stable even without that as- 
sumption. But it is slightly complicated to determine whether these are also sufficient or 
not. After a little lengthy calculations, we checked that if the conditions 



/4d < /4 min 



k 3 k n < kf 2 , 



k 2 

K d2 



k 2 
^dll 



^3^dl hikdi 



(D.10) 
(D.ll) 



are satisfied, Eqs. ( 1D.8|ID.9I) are also sufficient to ensure the stability of the above vacuum. 
The numerical solutions to the RG equations satisfy these conditions. 
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